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We consider a magnetic impurity in a spin liquid state of a magnetic system which is close to the 
quantum phase transition to the magnetically ordered state. There is similarity between this problem 
and the Kondo problem. We derive the impurity Green's function, consider renormalizations of the 
' magnetic moments of the impurity, calculate critical indexes for the magnetic susceptibilities and 

finally consider specific heat and magnetic interaction of two impurities. 
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I. INTRODUCTION 

<D 

. There is no need to explain the importance of the Kondo problem for condense matter physics. J± is relevant to the 
magnetic impurities in metals, heavy fermion compounds^!, tunneling phenomena in quantum dotscl, correlated lattice 
\ fermion systemsd, and many other physical systems. The problem of magnetic impurity in a two-dimensional (2D) 
insulating system with, long range antiferromagnetic order has also attracted great interest. This includes an impurity 
spin with an on-siteaa and sublattice symmetrical coupling, as well as an isolated ferromagnetic boncH. Because of 
Adler's relation for the impurity-spin-wave interaction these systems have no nontrivial infrared dynamics at zero 
temperature in spite of the gapless spectrum of Goldstone spin waveso. This makes the. impurity problem for the 
insulating state much simpler than the Kondo one. However it has been realized recentlya that in the case when the 
' 2D magnetic insulating system is close to the quantum critical point the infrared dynamics of the impurity is highly 
l— ', nontrivial and to a large extend is similar to that for the Kondo problem. The quantum critical impurity problem 
has been very recently addressed by Vojta, Buragohain, and Sachdevtli They demonstrated that there are nontrivial 
critical indexes for the impurity Green's function and for magnetic susceptibilities. In the present work we consider 
[ — i the same quantities calculating more accurately the critical indexes and in some cases the prefactors. In addition we 
I/"") ' consider specific heat and interaction of the impuritiesEa. 

o : 

O II. THE HAMILTONIAN 

o . 
o . 

To be specific we consider a magnetic impurity in two coupled Heisenberg planes. The two coupled Heisenberg 
planes is an ideal example of a 2D critical system which can be described by 0(3) nonlinear a- model. The two 
planes model is very well studied aij-ae«)±emperature both numerically&lij and analyticallylijO. Finite temperature 
I ■ properties are also well understoodt3'li3~E3. Hamiltonian of the system with impurity is of the form 

; 

H = H 2 + H imp , (1) 

H * = J H ( S W + s ? )s f ) + E Si w Si P0 . 

^ , <ij> i 

Here S s is spin 1/2 on the square lattice, the index i numerates cites, and the index n numerates planes. J is an 
antiferromagnetic coupling in the plane, and Jj_ is an antiferromagnetic coupling between the planes. rSpm of the 
impurity, s = 1/2, is coupled to one of the planes. It is known that in this system at J_l = 2.525 ± 0.002 Jl3 there is a 
quantum phase transition from quantum disordered state to the Neel state. In the present work we consider only the 
quantum disordered phase including the critical point. As we have already pointed out an interesting regime arises 
only close to the critical point and therefore we concentrate on the vicinity of this point. Note that in terms of the 
non-linear 0(3) er-modcl an effective Lagrangian of the system is L — d^d^ip — m 2 if 2 + ~f<ps, (p 2 = a 2 , where the 
parameters can be expressed in terms of parameters of the original Hamiltonian ([!]) . In our considerations we will use 
only the original Hamiltonian (Q). _ . 
Using bond operator representation^ 
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S^ 1 ' 2 ) = i(±t±tt-itt xt), (2) 

the two plane Hamiltonian H 2 from (|l|) can be rewritten in terms of the operators t; = {t itX ,t iyi t itZ ), and then 
diagonaliaf>d-|by a combination of the Fourier and Bogoliubov transformations with account of the hard core constraint, 
see Refs.EJO: 

U = ^e iqri (u q t q + « q t] 1 ), (3) 
q 

H 2 -> ^CJ q t q t q , 



where Mq,fq = y ^j- ± ^ are Bogoliubov parameters. The operator t q creates quasiparticle of the system. This 
quasiparticle has spin 1 and we call it spin wave or magnon. Near the critical point the spin-wave excitation energy is 

ujq « VA 2 + cV, (4) 

where A -C J is the spin- wave gap and c ~ 1.9 J is the spin- wave velocity. The function Aq in the vicinity of the 
critical point is q-independent: Aq w A w 2.4J. 

In the bond operator representation the impurity Hamiltonian is of the form 

#im P = |s (t + t f - it f x t ) -» |s (t + t f ) . (5) 

Here we have dropped the term to' x to- The matter is that we are interested in nontrivial long-range dynamics, but 
one can prove that all diagrams generated by the to ' x to term are infrared convergent, and therefore its contribution 

is less important. An alternative way is to redefine Hi mp as Hi mp — > | s ^So^ — So^ 2 ^, then the to^ x to term is 

canceled out exactly. Note that for an integer impurity spin the to ' x to term could be much more important because 
it can give a bound state of the spin wave with the impurity, and hence can lead to the full screening of the impurity. 
Using ([}]) we rewrite Hi mp in terms of quasiparticle operators t q . 



H lm p = ^"^(Uq + VqjtqV + h.C. « |^ ^ — + fc.C., (6) 
q V q V 1 

where a is the impurity Pauli matrix. Hereafter we set J = 1, so all energies are measured in units of J. 

III. THE IMPURITY GREEN'S FUNCTION AT ZERO TEMPERATURE 

Let us calculate the impurity self energy X shown in Fig. 1. 



FIG. 1. a) Single loop impurity self energy, b) Two loop noncrossing contribution to the self energy. The solid line corresponds 
to the impurity and the dashed line corresponds to the spin wave. 



In single loop approximation, see Fig. la, using eq. (q) we find 

1 d 2 q 



U!q(e - UJq) (2tt) 2 



(7) 



Close to the critical point w q = \J c 2 q 2 + A 2 , and hence simple integration gives for e = (i.e. at the position of the 
quasiparticle pole) 

E«(0) = -a 2 lnA. ( 8 ) 
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Here 



l07TC^ 

is dimensionless coupling constant, and A ~ 2 J is the ultraviolet cutoff. Similar calculation for the second order self 
energy shown in Fig. lb gives 

S«(0) = ~lnA. (10) 

As soon as a 2 / A <C 1 the second order self energy is small compared to the first order one, E^ 2 ) <C E^, and hence 
the perturbation theory is justified. However at a 2 / A > 1 the expansion does not converge and hence one has to sum 
all orders of perturbation theory. Exactly at the critical point, A = 0, the expansion diverges at arbitrary small a. 

The problem under consideration has a small parameter which is independent of the interaction. This is 1/N, 
where N=3 is number of components of the spin- wave excitation (O(N) a- model). In the leading in N approximation 
only the rainbow diagrams contribute to the impurity self energy. Summation of these diagrams leads to the usual 
noncrossing approximation (= self consistent Born approximation) for the impurity Green's function 

e- S(e) 



E(e) = a 2 J G{e-uj)dLo. (12) 

It is convenient to rescale the variables: e — > e/a 2 , ui — > Ld/a 2 , G — > a 2 G, E — > E/a 2 . In the new variables a 
disappears from the eqs. ( |ll| , |l2| ). Dependence on a remains only in the limits of cj-integration: E = J^ a Gdu>, where 

A Q = A/a 2 and A a = A/a 2 . Consider first the critical point, i.e. A a = 0. In this case eqs. ( |Tl|Jl^ ) can be solved 
analytically. The answer is 

G(e) = - ,2 v (i3) 



where eo is the impurity binding energy. Eq. (|13|) is valid if |eo — e| < 1. For illustration we present in Fig. 2 spectral 
function, — ilmG(e), obtained by direct numerical solution of eqs. ( |ll| , [i"2| ) at A Q = and K a — 1, 2, 5. 




FIG. 2. The impurity spectral function — — ImG(e) at the critical point. The curves obtained in the noncrossing approximation 
for different values of the ultraviolet cutoff: A a = 1, 2, 5. 
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Agreement with analytical solution ( |l3| ) is perfect. At the large A a = A/a 2 the impurity binding energy eo ~ —2.5. 
In the original variables it means that eo ~ —2.5a 2 , when a 2 <C 1. 

As one shall expect the Green's function ([D]) has no quasiparticle pole. However the pole appears away from the 
critical point when the spin-wave gap is nonzero . For illustration we present in Fig. 3 the spectral functions obtained 
by numerical solution of Eqs. (pl]Jl2[) for A a = 2, A Q = 0, and A Q = 2, A a — 0.05. 



1.5 
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FIG. 3. The impurity spectral function — — ImG(e) obtained in the noncrossing approximation. The dashed curve corresponds 
to the critical point (A a = 0), and the solid curve corresponds to the spin-wave gap A a — 0.05. In both cases the ultraviolet 
cutoff is A a — 2. 

It is clear that the quasiparticle peak absorbs spectral weight of the incoherent Green's function ( [13] ) from the area 
eo — e ^ A a . Therefore the quasiparticle residue Z oc v / 'A a (c.f. with Ref.EI). Slightly more detail analysis of the Eqs. 
(0121) shows that 



Zk, 0.8v/A a . 



(14) 



We have considered above the leading in N approximation. Let us estimate now 1/N correction which is due to the 
single loop contribution to the impurity-magnon vertex function shown in Fig. 4. 






FIG. 4. The impurity - spin wave vertex with account of single loop correction. The solid line corresponds to the impurity 
and the dashed line corresponds to the spin wave. 



Taking into account the algebraic relation for the Pauli matrices, cr^ovc^ = —a u , we find the vertex function given 
by Fig. 4 



T(e, A) = T bare (l - -1 J G(e - uj)G{e -u- \)dw 



(15) 



Here e and A are energies of incoming impurity and spin wave correspondingly, and Tb are is the bare vertex given by 
eq. (||) . Taking A ~ e — £q and using Green's function ([l3]) we find after integration in ( p"5| ) 



I\ = Th„ re [ 1 111 — - 

" ' 2N A 



(16) 
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This is the the first term of 1/N expansion, and keeping in mind scaling behavior we find 



r.ocA*, y«^4 (17) 



The impurity self energy is given by the diagram presented in Fig. 5. 

y ~~ ~~ n 



FIG. 5. The impurity self energy with account of vertex corrections. 

It differs from ( |l2| ) by the vertexes: the bare vertexes are replaced by the "exact" ones given by eq. JP7|). Effectively 
this introduces an additional factor ui 2y in the integrand in eq. (|l^). Assuming power behavior of the Green's function 

G(e) cx 1 —— (18) 

W (e - e) s V ' 

and performing integration we find 

S(e) cx (e - e) 2y ~ x+1 + const. (19) 
Substitution of this self energy into Dyson equation (^lj) gives the following condition of self consistency 

(e - e)~* oc (co - ef^" 1 . (20) 
Therefore the Green's function critical index with account of the leading 1/N correction is 

x = \ +y= i(l + i)«0.67 (21) 

Away from the critical point the quasiparticle pole appears in the Greens function. It absorbs spectral weight of 
the incoherent Green's function ([l8]) from the area eo — e < A a . Therefore the quasiparticle residue scales as 

ZcxA^, z = 1-x = \~V= 0.33. (22) 



IV. MAGNETIC MOMENTS OF THE IMPURITY, SUSCEPTIBILITIES 

Following Vojta, Buragohain, and Sachdev0 we consider two different types of the magnetic interaction. The first 
one is an interaction when the magnetic field h interacts only with the impurity 

H$ = -2sh. (23) 

The second case is homogeneous magnetic field 

H ( IP =-2sh-^2S i (n) h. (24) 

It is clear that in the first case the renormalized magnetic moment is proportional to the quasiparticle residue Z ', 
and hence, according to ( p2[ ) it scales as 

oc A z . (25) 

If we consider the system exactly at the_cxitical point, but at finite temperature, then the effective spin-wave gap is 
equal to the temperature, see e.g. Refs.E3'E3: 

A —> A T ps 0.962T. (26) 
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Together with (J^) this gives the following dependence of the impurity magnetic susceptibility on temperature. 

y « = ^! oc ^ s ,^ = ^_ (27) 



For homogeneous magnetic field the impurity magnetic moment is not renormalized because the interaction (24) is 
proportional to the total spin which is conserved. The magnetic moment is certainly redistributed over the volume of 
size r ~ 1/A T , but the value is conserved. For this reason the impurity susceptibility is given by usual the Curie law. 

Am = I 

Ximp rp' \ J 



This conclusion is not quite f^vial since the magnon cloud size r — > oo at T — > 0. Besides eq. Q2q ) does not agree with 
a conclusion from the papertJ. Therefore in the next section we also present a diagrammatic prove of our statement. 
It is rather technical section and a reader who is satisfied by the general arguments can skip it. 

V. CORRECTIONS TO THE IMPURITY MAGNETIC MOMENT 

In the present section only homogeneous magnetic field ( p4| ) is considered, therefore to simplify notations wa-omit 
the superscript (II). We restrict our consideration by single loop corrections, and follow the way used in Ref.E3 for 
calculation of the spin- wave magnetic moment. First we discuss a zero temperature case. Single loop self energy 
is given by eq. (g). This as a simple perturbation theory in a, but one can also consider this as a contribution to 
renormalization group equations with running coupling constant a and running infrared cutoff A. Anyway, correction 
to the quasiparticle residue due to (JsJ) is following 

<9£W a 2 

SZ °° 3A~ = - A < 29 ' 
The impurity magnetic moment is renormalized according to the diagrams shown in Fig. 6. 

y ' 1 , ' y 1 = ' 1 , 



FIG. 6. Zero temperature single loop corrections to the impurity magnetic moment. The cross denotes an external magnetic 
field attached to the corresponding line. 

First contribution comes from the correction to the quasiparticle residue: Sfie a = SZ. Straightforward calculation 
gives 5[iQb = — a 2 /(3A) and 5fi 6c = 4a 2 /(3A). Altogether this gives 

5fi 6 = S^ 6a + Sfi 6 b + Sfiec = 0. (30) 

So as one shall expect there is no renormalization of the magnetic moment. 

At finite temperature the relation (^0|) remains valid since the only thing we have to do is to replace A — > At- 
Howeverp-at finite temperature there are also additional diagrams which are due to the heat bath of the excited 
magnona 2 ^. First of all these are the two contributions to the self energy shown in Fig. 7 



a b 

FIG. 7. Temperature induced contributions to the impurity self energy. Dashed lines denote magnons from the heat bath. 



G 



Both contributions are proportional to the magnon mean occupation number 

1 



it _ r 



(31) 



but they are of the opposite sign and exactly cancel each other. So they do not influence the position of the quasiparticle 
pole. However these diagrams contribute equally to the quasiparticle residue 



n a du) a 2a 2 (°° dx 



5Z T = -2a 2 / =L = / —=—. (32) 

J w 2 A T A x 2 (e x -\) v ' 

This gives a correction to the inpurity magnetic moment. Other thermally induced corrections to the magnetic 
moment are given by diagrams presented in Fig. 8. 



-X- 



^1 Ills, 



abed 
FIG. 8. Thermal corrections to the impurity magnetic moment. Dashed lines denote magnons from the heat bath. The cross 
denotes an external magnetic field attached to the corresponding line. 

Straightforward calculation gives 5us a — J/igb ~ —37^7 dx/[x 2 (e x — 1)] , and 

Sfisc = Sfiad = Ji° dx/[x 2 (e x — 1)]. Total thermally induced correction to the impurity magnetic moment is 

equal to 

dfi T = SZ T + Sfi Sa + 5u sb + Su Sc + 5u sd = 0. (33) 

Together with eq. (|3p| ) this proves that the impurity magnetic moment is not renormalized. 

Thus, in spite of the cloud of virtual and thermal magnons, effectively the impurity in the external magnetic field 
can be described as a spin 1/2 system with unrenormalized magnetic moment. Hence we immediately come to eq. 
(p8|). We would like to stress that the arguments related to the magnetic moment guarantee only singular in T part 
of the susceptibility, therefore instead of ( pgj ) it is more correct to write Ximp — ^/T + const. 



VI. SPECIFIC HEAT RELATED TO THE IMPURITY 



Binding energy of the impurity has been calculated in section III. Exactly at the critical point, at zero temperature, 
and at small coupling constant, a 2 ^ 1, the binding energy is e ~ —2.5a 2 . A finite spin-wave gap A pushes the 
position of the quasiparticle pole up, see Fig. 3. For a small gap the dependence of eo on the gap is linear, and the 
coefficient of the proportionality is approximately 2: 

e ~ a 2 (-2.5 + 2A). (34) 

The origin of A in this equation is not important: whether A is nonzero because the system is away from the critical 
point, or the system is at the critical point, but A is nonzero because of temperature. In the later case we can use 
eq. ( |26|) for the gap and hence the impurity specific heat is 

Cimp = ~ 2a 2 . (35) 

Note that this is highly unusual result because Ci mp ^ at T — + 0. For comparison: the bulk specific-beat of a 2D 
antiferromagnet at the quantum critical point is quadratic in temperature, C bu ik oc T 2 , see e.g. Refs£3'E3. From (|3q) 
one concludes that the impurity entropy is Si mp = J Ci mp dT /T cx InT. Strictly speaking this is nonsense because it 
gives Si mp (T — 0) = —00. This probably indicates that there is a small nonzero critical index £ in the specific heat 
dependence: Ci mp cx T^. Equation ( |34"| ) has been derived in the noncrossing approximation. One can check that single 
loop vertex corrections considered at the end of the section III do not change this equation. It probably means that 
nonzero £ can appear only due to the two loop corrections, £ cx 1/N 2 . There is no doubt that this is a very interesting 
problem which deserves further analytical analysis and which can be also studied in numerical simulations. 
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VII. STATIC INTERACTION BETWEEN TWO DISTANT IMPURITIES 



Interaction between two impurities is given by diagram shown in Fig. 9. 




FIG. 9. Spin-wave exchange between two impurities. 



Since the impurities are localized, we have to integrate over all possible momenta transfer q. Away from the critical 
point there are two distinct regimes: c/r < A and c/r > A, where r is distance between the impurities and c ~ 1.9 J 
is the spin-wave velocity. In the first regime the interaction drops down exponentially with distance and it is not an 
interesting case. We consider only c/r > A case which is also relevant to the critical point. We also assume that the 
coupling constant is strong enough: a 2 /A > 1. The impurity-magnon interaction is given by the Hamiltonian (j^). 
According to the consideration in section III it is renormalized as Hi mp — ■» Z q T(q), where Z q is quasiparticle residue 



([14), and T(g) is the vertex (17). Therefore the interaction corresponding to the diagram Fig. 9 is of the form 



rr/ v f Z IA^ d\ _ _ Af f [Z q T q /T bare ] 2 cfq 

V(r) = / -z-z — — - = -a x ■ o x — \ — — — — -, (36) 

where a\ and 02 are the impurities Pauli matrixes. One power of w q in the denominator appears because of the 
spin- wave propagator and another power is due to the bare interaction (^j). With account of (^), and (|l7|) this 

gives 

V(r) ~ -ox ■ ox a 2 -^ [ q 2 ^^e l ^d\ ~ (37) 



This is a long range ferromagnetic interaction. An interesting fact is that the interaction is independent of the bare 
coupling constant a. Another interesting fact is that there is no renormalization of the power in the 1/r dependence. 
However we stress that these facts have been proven only in the one loop approximation (one loop above the leading 
noncrossing approximation). 



VIII. CONCLUSION 



The dynamics of the magnetic impurity in a 2D antiferromagnet close to the critical point is highly nontrivial. To 
some extend it is similar to the dynamics in the Kondo problem. In the present paper we have considered spin 1/2 
impurity. We have calculated indexes for critical behavior of the Green's function, vertexes, magnetic moments and 
magnetic susceptibilities. We have also considered the impurity specific heat and long range interaction between two 
impurities which is due to the spin-wave exchange. 
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